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1.  Introduction. 


regressor,  a,  and  3 unknown  regression  coefficients  and  e a normal  ran- 

2 

dom  variable,  with  zero  mean  and  variance  a . For  a fixed  value,  1), 
x(H)  = (11-  oO/p  is  called  the  inverse  (linear)  regression  parameter. 

Given  n points  (x^,Y^),..., (xn,Yn)  one  can  determine  exact  y- level 
confidence  intervals  to  the  inverse  regression,  on  the  basis  of  the  least- 


the  American  Statistician,  Zerbe  (1978)  applied  Fiedler' s Theorem  to  ob- 
tain confidence  intervals  for  the  ratio  of  arbitrary  linear  combinations 


of  the  parameters  of  the  general  linear  model.  Zerbe' s generalization 
yields,  as  a special  case,  confidence  intervals  for  the  reltive  potency 
of  dilution  bioassays.  These  confidence  limits  sure  obtained  from  the  con- 


fidence limits  of  the  difference  of  two  inverse  linear  regression  para- 
meters, which  correspond  to  parallel  regression  lines.  This  however,  was 


of  non  - parallel  linear  regressions  one  cannot  obtain  exact  confidence 
intervals  as  a special  case  of  Zerbe' s method.  Formulae  available  in  the 
literature  are  generally  based  on  asymptotic  approximations  (see  Armitage 
(1971),  pp.  446-447).  In  the  present  paper  we  develop  a method  for  de- 
termining exactly  confidence  intervals  for  any  linear  combination  of  the 


corresponding  inverse  linear  regression  parameters  of  two  non- parallel 
lines*  The  method  is  based  on  a simultaneous  application  of  the  basic 


statistic  used  in  Fieller's  Theorem  to  the  bivariate  case.  We  obtain 
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a joint  confidence  set,  at  level  y , for  the  two  inverse  linear  regres- 
sion parameters.  The  set  is  closed  and  convex,  with  continuous  boundary, 
similar  to  an.  ellipse  which  has  been  asymmetrically  extended.  This  simul- 
taneous confidence  region  yields  simultaneous  confidence  intervals  for 
any  linear  combination  of  the  parameters,  by  considering  the  proper  tan- 
gential lines.  The  method  is  developed  in  Section  2 in  a general  manner 
and  reduced  by  proper  transformations  to  a canonical  form.  Fortran  sub- 
routine programs  for  determining  the  upper  and  lower  confidence  limits 
of  the  difference  axe  provided  in  Appendix  I.  In  Section  3 we  compare  the 
exact  simultaneous  confidence  limits  of  the  difference  to  the  asymptotic 
limits.  The  comparison  is  done  numerically  by  simulating  samples  of 
size  n = 50  and  n = 100;  determining  the  least- squares  estimators  of 
the  regression  lines,  the  inverse  regression  estimates  and  the  confidence 
limits  for  the  difference  of  the  inverse  regression  parameters.  It  is 
shown  that  the  asymptotic  approximations  yield  intervals  which  are  too 
short  compared  to  the  exact  method,  even  in  samples  of  size  n = 100. 

This  illustrates  that  the  commonly  applied  asymptotic  formulae  lead  to 
intervals  with  smaller  coverage  probability  than  the  nominal  confidence 
level. 

There  is  in  the  literature  a considerable  disagreement  on  the  ques- 
tion whether  the  notion  of  relative-potency  is  meaningful  in  non-dilution 
assays.  Cornfield  (1964),  in  his  famous  paper  on  the  role  of  parallelism 
in  comparative  bioassays,  provided  a long  discussion  of  the  theoretical 
aspects  of  this  question.  Cornfield  developed  in  his  paper  a Bayesian 
procedure  for  determining  (posterior)  confidence  intervals  of  the  log- 
relative-  potency,  as  functions  of  the  inverse  regression  of  the  test 

i 

i 
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preparation,  xt(T)).  We  apply  in  Section  4 the  method  developed  in  the 
present  paper  to  obtain  classical  confidence  intervals  of  log-relative- 
potency  of  comparative  (non- dilution)  photodynamic  bioassays.  These 
assays  were  performed  by  Epstein  et  al  (1965)  to  test  the  toxicity  of 
benzen- soluble  organic  extracts  from  air  pollution  samples.  The  standard 
preparation  consisted  of  factory  made  benzo  (a)  pyrene  (BaP).  Some  de- 
scription of  the  assays  and  the  nature  of  the  biological  response  is  given 
in  Section  4.  We  illustrate  the  method  developed  here  on  actual  sample 
data  from  Boston,  Mass,  and  Chattanooga,  Tenn.  For  an  extensive  statis- 
tical analysis  and  new  indices  based  on  photodynamic  bioassays  see 
Bialik,  Epstein  and  Zacks  (1978)  and  Bialik  (1978). 

In  conclusion,  the  method  developed  here  of  determining  exact  confi- 
dence intervals  is  shown  to  be  easily  applicable.  Although  it  requires 
computer  analysis,  the  program  is  simple  and  the  results  can  be  obtained 
in  a matter  of  seconds  on  any  time- sharing  equipment.  Programmable  cal- 
culators can  be  vised  too  for  determining  these  confidence  limits. 

2.  Derivations. 

Consider  two  regression  lines: 

Y = a + 8 x + e (the  "standard"  line) 

and  s s 

Y * + Ptx  + e (the  "test"  line) 

where  ag,  at  and  are  the  regression  coefficients  and  e is  a 

random  variable  having  a normal  distribution,  N(0,o2).  It  is  assumed 

2 

that  the  variance,  0 , around  the  two  lines  is  the  same.  n indepen- 
dent  trials  are  conducted  to  estimate  the  coefficients  of  the  "standard" 


1 
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line, and  independent  trials  are  conducted  to  estimate  the  coefficients 

of  the  "test"  lines.  Let  {xe  } be  the  points  at  which  the 

s,  ± s } ng 

"standard"  trials  are  performed  and  {x.  -,..., x.  ] those  for  the  "test" 

1 t,  n-t 

trials.  We  denote  by  xo  and  x.  the  averages  of  x . (i=l, ...,n  ) 

s « S|1  ' s 

ns  _ 

and  x.  . (l=l,...,n  ),  respectively.  Let  -SD  = S (x  , - x y and 
“y1  s i=l  s,  I s 

define  SD^  similarly.  It  is  assumed  that  the  "standard"  and  the  "test" 

trials  yield  mutually  independent  random  variables  Y . and  Y. 

s,  i t,  J 

(i  = j = 1,  ...jn^). 

Denote  by  a , b , a.  and  b.  be  the  least- squares  estimators  of 

S S u o 

a2 

as*  at  811(1  respectively.  Let  o be  the  pooled- variance 

around  the  least- squares  regression  lines  (see,  Graybill  (1977)  for  the 
theory  and  formulae  of  least- squares  estimation).  Fix  and  let  x 

s 

and  Xj.  denote  the  estimators  of  the  inverse  regression  parameters 
§s  = (^“asVPs  and  = (T)“  c»t)/Pt,  in  which  the  least- squares  esti- 
mators a , b ; a and  b are  substituted,  respectively.  We  define 
the  pivotal  variables 


u = 


(H  - *s) 


5b  . b (S  - 5), 
3s  s s s 3S  * 


(2.1) 


- s’  - Vt  ’ bt(it  - V • 


It  follows  immediately  from  the  least- squares  theory  that  u and  u. 

s t 

are  independent  random  variables,  normally  distributed  with  zero  means 
and  variances 

k2 


2 /,  (x  - 5 V 

Var{u  ) = a2(—  + ■ s_  s 
aJ  In  SD 


2/i  GL- 5t)‘ 

V- 


(2.2) 


f 
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It  follows  that, 

(2.3) 


.2/* 
n b (x 
s s s 


V 


- V2 


1 + 


2 2r_  i 

o X [2] 


SD 


2 

where  x Cv3  designates  a chi- squared  random  variable  with  v degrees  of 

a2 


freedom.  'Hie  pooled-variance  estimator  a is  independent  of  the  random 
variable  specified  in  (2.3)  having  a distribution  like  that  of  oZ%[\>  ]/v  , 


where  vp  = ns  + nt  “ 4.  ] denote  the  yfractile  of  a ( 


P ' P' 


cen- 


tral) F- distribution  with  2 and  degrees  of  freedom.  It  follows 


that 


(2.4) 


,2/a  _ \2 

VV*.’ 


1 + 


h (^ -§  )‘ 
s s s 


ifrt'VV*  2.2 

nt(xt-?t)2  * FyC2,Vp] 


SD 


1+- 


SD, 


holds  with  probability  y for  all  (a,, Pg,at, Pt,o).  Inequality  (2.4) 

specifies  a y- level  simultaneous  confidence  region  for  (5  ,5.).  Define, 

s t» 


(2.5) 


sf  (5.  - 5J  » 

S 


SDt  " V ' 


and 


(2.6) 


1 


l 


□ □ 


(2.7) 


n s 9 
8 


a** 


then  inequality  (2.4)  is  reduced  to  the  canonical  form 


(2.8) 


Bs**s  - C8>  , Bt<*t  ’ Ct>  , 


i + c: 


i + c; 


2 *2 

where  R * 2a  ?^[2,v^]/ns.  Since  both  summands  on  the  LHS  of  (2.8)  are 
non- negative,  it  follows  that  a necessary  condition  for  (Q  ,Q.)  to  satis- 
fy  (2.8)  is  that  s(ss  Cg2),  "here  (i  = 1,2)  are  the  two 

real  solutions  of  the  quadratic  equation 


(2.9) 


- <s>  - * 0 . 


These  roots,  if  exist,  are  given  by 


(2.10)  C 


s s , -»i  R r-2/.  ^ -2%  „2,i 

b2  _ r2  + ^ b2  _ r2[Bs(1  + zs)  " R ) • 


2 1 

notice  that  R = 0(— ) and  therefore,  if  n is  sufficiently  large, 

2 2 s , 99 

B > R with  high  probability.  Indeed,  SD  /n  = 0(l)  and  b -*  Br  > 0, 
8 S'  s s s r 

with  probability  one.  If  the  real  roots  (2.10)  do  not  exist,  (2.8)  does 
not  hold  and  one  cannot  obtain  exact  confidence  limits.  The  probability 
of  such  an  event  is,  however,  smaller  than  1-y.  For  every  Cg  In  the 
interval  (C^C^),  define 

3 3 


(2.11) 


WM  * R - BX  - V 7(1  ♦ c2). 


i 
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It  follows  from  (2.8)  that  the  corresponding  value  of  C*  for  a point 
vithin  the  confidence  region,  is  between  the  two  limits 


.(i) 


Vt 


(2.12) 


Ct  <Cs>  T*.  vc,.  j * c-1)i[wssK  • a * *»> 

« R s s 


In  Figure  1 we  present  a computer  graphing  of  the  boundary  (2.12), 
for  the  case  of  n^=  ns  = nt  = 20,  zg  = .4385,  ^ = -.390,  Bs  = I.9l8, 

t 7 and  R _ .049917.  The  confidence  region  prescribed  by  the 

boundaries (2. 10)  - (2.12)  is  convex  and  has  a smooth  boundary.  It  is 
specified  by  a fourth  degree  curve,  which  looks  like  an  asymmetrically 
extended  ellipse.  Simultaneous  confidence  intervals  for  any  linear  com- 
ination  u)  _ X^s  + XgC^,  with  specified  values  of  and  Xg,  can 
be  obtained  by  considering  the  two  tangential  lines  having  slopes 

intercePts  of  thes*  lines  provide  upper  and  lower  confidence 
limits  for  o,/^.  In  particular,  for  the  difference  6 = Ct  - C , we 

consider  the  two  tangential  lines  with  slope  1 and  intercepts  \ and 

which  are  the  upper  and  lower  confidence  limits  for  6.  For’  every 

* iQ  the  interval  (C<1,,C<2))  consider  the  functions. 


(2.13) 
and 

(2.14) 


U(x 


) _ Bt2t  * tTR^«s)Bt(l^-fg(xj2J]* 

I2  T ~ 


L(x) 


*t  - V*2s> 

Bt2t~ 


Bt  “ TR^zs^ 


1 


FIGURE  1.  Computer  Graphing  Of  The  Boundary  (2.12)  For  n = n = 20, 

S T» 
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These  two  functions  consist  of  the  upper  and  lower  boundary  points  of  the 
confidence  region  defined  by  (2.12).  The  function  U(x)  is  concave,  having 
a maximum  at  x = zg.  Let  xu  be  the  (unique)  point  at  which  U(x)  = 1, 
and  ^ the  point  at  which  ~ L(x)  = 1.  It  is  easy  to  verify  that 


(2.15) 


4 - 


Vt 


dx  | / _j2  / A \ \ 2 

• \S\  “ YR  X*  Zs 

i[YR(x^S)Bt(l  + - YI(x^s)]_i 


B*(l  + zl)  - 2Y  (x,z  ) 

— ^ ^ + [Y  (xjz  )B*(l  + z2) 

B2  - Y (x;z  R s t t 


where 


(2.16) 


" VXiZs^  / ^Bt  " VX'Zs^  ( 3x  VX'Zs^ 


a 2B^(z  - x) 

dx  VX'Zs^  = , Z\2 

(1  + x ) 


The  derivative  of  L(x)  can  be  expressed  similarly,  with  minus  signs  be- 
fore the  second  and  third  terms  in  the  curly  brackets  of  (2.15).  The 
solution  of  the  equations  U(x)  = 1 and  L(x)  = 1 is  performed 
numerically.  The  upper  and  lower  confidence  limits  for  6 are  given  by 


(2.17) 


6u  - “<V  - *u  ’ 


h ■ u(*l>  - *1 


In  Appendix  I we  provide  the  subroutine  functions  which  determine  numer- 
ically the  values  of  5 and  &r  . 

U L 
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We  conclude  the  present  section  with  some  comments  about  the  asymp- 
totic (large  sample)  confidence  intervals.  It  can  be  shown  (see  Zacks 

(1971),  pp.  24l)  that  the  asymptotic  distribution  of  z and  z.  axe 

s t 

normal  with  means  £s  and  , respectively;  and  asymptotic  variances 

(2.18) 


and 


AV(2S)  - -f-  (1  + C.)  , 

s Bn  s 

s s 


AV{£  ) * (1  + cp  . 

Accordingly,  the  asymptotic  confidence  limits  of  6 = Q.  - £ , at  nominal 

u S 

level  Y,  are  given  by 

T .2  2"l* 

* - 1 + z 1 + zf 

(2.19)  Si  = (zt  - zs)  + (-1)1  U + — _ , i =1,2 

-5-  V ns  B B. 

2 L S t 

where  Uj.  is  the  (i^)th  fractile  of  the  standard  normal  distribution. 

~ 

In  the  following  section  we  compare  the  behavior  of  the  asymptotic  limits 


with  those  of  the  exact  ones. 


Simulation  Comparison  of  the  Exact  and  Asymptotic  Confidence 
Intervals  for  Differences. 

In  the  present  section  we  illustrate  the  exact  and  the  asymptotic  con- 
fidence intervals  by  performing  simulation  runs  and  computing  the  confidence 
limits  for  the  simulated  values.  For  the  sake  of  simplicity  we  consider 
two- point  designs,  in  which  both  the  "standard”  and  the  "test"  are  performed 
at  the  points  z = -1  and  z = 1.  Altogether  there  are  n = 2m  observa- 
tions around  each  regression  line,  m observations  at  z =-l  and  m at 
z *»  1.  When  there  are  only  two  regression  points  the  least- squares  re- 
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gression  line  passes  through  the  sample  means  corresponding  to  these  points. 


Thus,  we  have  generated  in  each  ca'e  two  independent  standard  normal  vari- 


ates U and  simulated  the  sample  means  by 


*-i  ■ “-i  + "5=  « 


^ + ^ ■ 


The  slope  of  the  regression  line  is  then  b = (y^  - Y ^)/2  and  its  inter- 
cept is  a = (Y  + Y ^)/2. 


After  computing  these  estimates  for  the  "standard"  and  the  "test" 


regression  lines  we  determined  z = (T]~  a )/b  , z.  = (11-a  )/b  and  the 

S S S b b *b 


exact  confidence  intervals,  at  level  y = .95#  according  to  the  subroutine 


functions  of  Appendix  I.  In  our  computations  we  assumed  that  a is 


2 2 2 

known  and  applied  the  formula  R = a ^[2]/n.  Asymptotic  confidence 


limits  were  computed  in  each  run  according  to  formula  (2.19)  using  the 


time  value  of  a instead  of  the  estimator  5.  In  Tables  1 and  2 we  pre- 


sent the  results  of  50  such  simulation  runs..  Table  2 presents  cases 


with  sample  sizes  twice  as  large  (n  = 100)  as  those  of  Table  1 (n  = 50). 


We  see  that  the  asymptotic  limits  yield  intervals  which  are  too  small. 


even  when  n = 100.  Thus,  the  asymptotic  formulae  yield  intervals  whose 


coverage  probabilities  are  smaller  than  the  nominal  ones.  We  remark  that 


the  computation  of  50  runs  took  about  six  seconds  on  a time- sharing 


(Honeywell  GE-U30)  computer.  The  determination  of  the  exact  confidence 


limits  in  each  case  is  a matter  of  a split  of  a second. 


l 


TABLE  1.  .9 5- Confidence  Limits  for  5j  n = 50;  a =3.5, 

s 

at=3.5,  = 2. 5 1 P^.=l»75 1 = 0,  0 = .5  . 


Point  Estimates 

Exact 

Limits 

A • 

z 

A 

z. 

A 

5 

S_ 

A 

5 

s 

t 

L 

u 

0. 141 

0.  10£ 

-.  039 

-.  199 

0.  1£3 

0.  £75 

0.  £09 

-.  066 

-.££8 

0.  099 

0.  185 

0.  146 

-.  039 

-.  £Q£ 

0.  1£7 

0.  178 

0.  088 

-.  090 

-.£45 

0.  067 

0.  179 

0.  143 

-.  036 

-.  191 

0.  1£0 

0 . £4  0 

0.  £15 

-.  0£5 

-.  187 

0.  140 

0.  £55 

0.  1£1 

-.  134 

-.£86 

0.  019 

0.  £14 

0.  £04 

-.  009 

-.  17£ 

0.  157 

0.  ££0 

0.  186 

-.  034 

-.  191 

0.  1£7 

0.  £14 

0.  £13 

-.  00£ 

-.  156 

0.  156 

0.  150 

0.  148 

-.  00£ 

-.  153 

0.  150 

0.  £03 

0.  145 

-.  058 

-.£14 

0.  101 

0.  167 

0.  130 

-.  037 

-.  188 

0.  116 

0.  £15 

0.  169 

-.  046 

-.  199 

0.110 

0.  194 

0.  188 

005 

-.  165 

0.  157 

0.  £18 

0.114 

-.103 

-.£55 

0.  049 

0.  £41 

0.  079 

— .1 6£ 

-.338 

0.  015 

0.  191 

0.  094 

-.  097 

- . £5£ 

0.  058 

0.  £48 

0.  069 

-.  178 

-.335 

-.  0££ 

0.  188 

0.  186 

-.  00£ 

-.  171 

0.  17£ 

0 . £64 

0.  178 

-.  086 

-.  £49 

0.  080 

0.  153 

0.  £03 

0.  051 

-.111 

0.  £18 

0 . £ 1 6 

0.  134 

-.  083 

-.£36 

0.  07£ 

0.  154 

0.103 

-.  051 

-.  £1£ 

0.114 

0.  15£ 

0.  06£ 

-.  090 

-.£40 

0.  060 

0 . £ 1 £ 

0.  £41 

0.  0£9 

-.  1£5 

0.  187 

0.  183 

0.  060 

-.  1£3 

-.£87 

0.  041 

0.  177 

0.  169 

-.  008 

-.  173 

0.  161 

0.  1S£ 

0.  1£5 

-.  057 

- . £ 1 £ 

0.  099 

0.  160 

0.  093 

-.  068 

-.££5 

0.  091 

0.  193 

0.  161 

— . 03£ 

-.  181 

0.  119 

0.  183 

0.  073 

-.110 

-.£76 

0.  057 

0.  £39 

0.  133 

-.106 

-.£65 

0.  055 

0.  £08 

0.  £04 

-.  004 

-.  161 

0.  15S 

0.  £33 

0.  £3£ 

-.  001 

-.  1 7£ 

0.  177 

0.  ££8 

0.  173 

-.  055 

-.£15 

0.108 

0.  ££4 

0.  070 

-.  154 

-.313 

0.  004 

0.  £35 

0.  199 

-.  037 

-.£04 

0.  135 

0.  £35 

0.  071 

-.  165 

-.  3£3 

-.  006 

0.149 

0.  197 

0.  048 

-.  1£1 

0.  £££ 

0.  £1 7 

0.  196 

-.  0£1 

-.  185 

0.  148 

0.  £47 

0.  1££ 

-.  1£4 

-.£91 

0.  045 

0.  ££0 

0.  061 

-.  160 

-.318 

-.  000 

0.  ££4 

0.  1 10 

-.114 

-.£73 

0.  047 

0.  £04 

0.304 

0.100 

-.  068 

0.  £74 

0.  184 

0.  153 

-.  030 

-.  191 

0.  133 

0.  £56 

0.  181 

-.  076 

-.£35 

0.  087 

0.  £66 

0.  181 

-.  085 

-.£41 

0.  07£ 

0.  £69 

0.  145 

-.  1£4 

-.£87 

0.  040 

0.  ££3 

0.£17 

006 

-.178 

0.  171 
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TABLE  2.  .9 5- Confidence  Limits  for  6}  n = 100,  ag  » 3.5, 

at  = 3*^>  Ps  = = 1*^5,  11  * 4.0,  o » .5  . 


Point  Estimates 

Exact 

Limits 

2 

s 

* 

zt 

A 

5 

a 

bl 

6 

u 

0. 168 

0.121 

-.048 

-.  134 

0.  040 

0.241 

0.  179 

-.  062 

-.  150 

0.  026 

0.  192 

0.  144 

-.  047 

-.  135 

0.  041 

0.  188 

0.  113 

-.  075 

-.  160 

0.  011 

0.  189 

0.  143 

-.  046 

-.  131 

0.  040 

0.221 

0.  183 

-.  039 

-.  126 

0.  049 

0.230 

0.  131 

-.100 

-.  184 

-.  015 

0.207 

0.  177 

-.  031 

-.118 

0.  058 

0.211 

0.  166 

-.  044 

-.  131 

0.  042 

0.208 

0.  182 

-.  026 

-.111 

0.  060 

0.  173 

0.  146 

-.  027 

-.110 

0.  058 

0.201 

0.  144 

-.  057 

-.  143 

0.  029 

0.  181 

0.  136 

-.  046 

-.  129 

0.  039 

0.208 

0.  157 

-.  051 

-.  135 

0.  035 

0.  197 

0.  168 

-.  029 

-.115 

0.  053 

0.210 

0.  127 

-.  083 

-.  167 

0.  002 

0.222 

0.110 

-.112 

-.203 

-.  020 

0.  195 

0.116 

-.  079 

-.  164 

0.  006 

0.225 

0.102 

-.  124 

-.209 

-.  033 

0.  193 

0.  165 

-.  028 

-.118 

0.  062 

0.234 

0.  162 

-.  072 

-.  160 

0.  016 

0.  174 

0.  175 

0.  001 

-.  086 

0.  090 

0.209 

0.  138 

-.  071 

-.  156 

0.  014 

0.  174 

0.  122 

-.  053 

-.  140 

0.  035 

0.  174 

0.  098 

-.  076 

-.  159 

0.  008 

0.206 

0.  198 

-.  009 

-.  094 

0.  077 

0.  191 

0.  099 

-.  092 

-.  179 

-.  004 

0.  188 

0.  157 

-.  031 

-.119 

0.  059 

0.  190 

0.  133 

-.  057 

-.  143 

0.  028 

0.  178 

0.  116 

-.  062 

-.  148 

0.  024 

0.  196 

0.  153 

-.  043 

-.  126 

0.  041 

0.  190 

0.106 

-.  085 

-.  173 

0.  004 

0.222 

0.  138 

-.  084 

-.  170 

0.  003 

0.204 

0.  176 

-.  028 

-.114 

0.  059 

0.218 

0.  189 

-.  029 

-.119 

0.  062 

0.215 

0.  159 

056 

-.  143 

0.  032 

0.213 

0.103 

-.110 

-.  196 

-.  024 

0.219 

0.  172 

-.  047 

-.  136 

0.  043 

0.219 

0.  103 

-.116 

-.202 

-.  030 

0.  172 

0.  171 

-.  001 

-.  090 

0.  089 

0.209 

0.  171 

-.  038 

-.  126 

0.  051 

0.226 

0.  132 

-.  094 

-.  182 

-.  004 

0.211 

0.  098 

-.113 

-.  199 

-.  026 

0.213 

0.  125 

-.  088 

-.  174 

-.  001 

0.202 

0.230 

0.  028 

-.  061 

0.118 

0.  191 

0.  149 

-.  043 

-.  130 

0.  045 

0.231 

0.  163 

-.  068 

-.  154 

0.  020 

0.237 

0.  164 

-.  072 

-.  157 

0.  014 

0.238 

0.  144 

-.  094 

-.  181 

-.  006 

0.213 

0.  181 

-.031 

-.  121 

0.060 

w*  ■ . ! • \ 

“ • •■'*»  ■ 

• • • ,*r 

— ■ - — - 

Asy.  Limits 


0.  008 


0.  025 


0.  022 

0.  018 


0.  044 


0.  072 


0.  018 


0.  060 


0.  041 
0.  011 
0.  007 
0.  024 


0.  042 
0.  044 
0.  014 


0.  025 


0.  033 


0.  100 
0.  027 
0.  002 
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4.  Confidence  Intervals  for  the  Relative- Potency  in  Comparative  Bloassays. 

In  the  present  section  we  illustrate  the  application  of  our  method  to 
the  determination  of  exact  confidence  limits  of  the  log- relative- potency  in 
comparative  bioassays.  The  experiments  described  here  were  performed  in 
1963  and  1964  (see  Epstein  et  al  (1965)).  They  consisted  of  a series  of 
photodynamic  bioassays  in  which  air- pollution  samples  from  various  sites 
in  the  U.S.A.  were  chemically  analyzed.  The  benzo- soluble  organic  extracts 
were  then  diluted  in  aceton  and  the  substances  applied  on  a Paramecium 
Cauda turn  under  ultraviolate  radiation.  The  measurement  of  response  was 
the  time  required  (in  minutes)  to  kill  (or  deactivate)  9056  of  30  cells  in 
the  suspension.  This  observed  random  variable  is  called  the  LT90.  Three 
dosages  of  the  organic  extracts  were  applied  (10-^,  10~  and  lcf  ^ [g/ml]). 
It  was  empirically  shown  that  Y = dm.  (LT90)  is  approximately  normally 
distributed,  and  that  the  expectation  of  Y is  related  to  the  log-dose 
linearly.  Each  assay  was  repeated  m = 4 times  independently.  A standard 
assay  was  also  performed  with  similar  doses  of  factory  made  Benzene- a- Peryne 
(BaP).  Also  in  the  standard  preparations  Y = tm.  (LT90)  was  linearly  re- 
lated to  the  log-dose  of  the  BaP.  The  slopes  of  the  regression  lines 
were,  however,  significantly  different.  The  log-relative-potency  at  level 
H is  defined  as  the  difference  6 between  the  inverse  regression  para- 
meters §s  and  §t.  In  Table  3 we  provide  the  LT90  values  and  the  least- 
squares  estimates  of  the  regression  parameters.  We  include  the  "test" 
preparations  from  assays  performed  with  the  1963- air- pollution  data  from 
Boston,  Mass,  and  Chattanooga,  Term.,  and  corrsponding  "standard"  prepares 
tions.  In  Table  4 we  provide  the  exact  .95“ confidence  limits  for  the 
log- relative- potency,  corresponding  to  LT90  = 30(2)  50.  These  confidence 


15- 


f 


limits  vrere  computed  according  to  the  subroutine  functions  of  Appendix  I, 
after  reducing  the  statistics  to  the  canonical  ones.  We  see  in  Table  4 
that  corresponding  to  LT90  of  40  minutes,  the  standard  preparation  of  BaP 
is  at  least  a hundred  thousands  times  more  potent  that  the  organic  extracts 
from  Boston  and  at  least  two  million  times  more  potent  than  those  from 
Chattanooga.  The  same  method  can  be  used  to  compute  the  confidence  limits 
of  the  relative  potency  of  the  organic  extracts  from  Chattanooga  compared 
to  Boston.  These  confidence  limits  are  presented  in  Table  5. 


TABLE  3.  The  LT90  Values  and  the  0n  (LT90)  to  log-dose 
Regression  Estimates,  for  the  Photodynamic  Bio- 
assays. 


Boston,  Mass. 

Chattanooga,  Tn. 

BaP 

dose  [g/ml] 

BMiia 

ESKEHESI 

10" 4 io“5  10“ 6 

LT90  (min) 

10.91  25.86  85.62 

14.93  37.42  85.23 

5.97  7.4l  9.37 

14.38  24.59  83.21 

15.23  38.21  82.41 

5.68  6.83  8.97 

14.22  28.53  86.50 

15.68  38.78  83.16 

6.22  7.38  9.64 

.14.68  28.81  87.04 

15.79  39.01  84.02 

7.04  7.64  9.80 

a 

3.4468  + .06324 

3.6029  + .03162 

2.0199  + .03162 

b 

.9251  + .07746 

.8463  + .03873 

.2096  + .03873 

o2 

.012 

.003 

*)  - 

*)  a is  the  pooled  estimate  of  the  variance  around  the  "standard" 
and  the  "test"  regression  lines. 
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TABLE  4.  .95-Confidence  Intervals  of  the  Log-Relative  Potency 
of  the  Standard  (BaP)  Against  the  Organic  Extracts. 


LT90 

Boston 

Chattanooga 

Lower 

Upper 

Lower 

Upper 

30.0000 

4.3784 

13.5967 

5.5451 

8.9085 

32.0000 

4.5138 

14.1597 

5.7176 

9.2310 

34.0000 

4.6398 

14.6871 

5.8798 

9-5413 

36.0000 

4.7588 

15.1839 

6.0328 

9.8334 

38.0000 

4.8720 

15.6530 

6.1778 

10.1086 

40.0000 

4.9799 

16.0960 

6.3157 

10.3688 

42.0000 

5.0804 

16.5135 

6.4473 

10.6153 

44.0000 

5.1776 

16.9221 

6.5735 

10.8492 

46.0000 

5.2708 

17.3150 

6.6906 

11.0749 

48.0000 

5.3586 

17.6859 

6.8028 

11.2952 

50.0000 

5.4442 

18.0355 

6.9111 

11. 5060 

TABLE  5.  .95-Confidence  Limits  to  the  Relative  Potency  of  the 
Organic  Extracts  from  Chattanooga  Compared  to  Boston. 


LT90 


30 

36 

42 

50 


Lower  Limit 


.471 

.494 

.507 

.513 


Upper  Limit 

.883 

.918 

.963 

1.032 
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Appendix  I 

We  present  here  subroutine  functions,  in  FORTRAN,  for  the  computation 
of  the  upper  and  lower  confidence  limits  of  the  difference  8 = Q.  - Q . 

w S 

2 2 a a 

The  input  variables  are:  A «-  B , B«-B,  x«-z,  Y*-z.  and  R ♦-  R . 

Swob 

*2 

The  confidence  level,  sample  size,  and  the  value  of  a are  all  deter- 
mining the  value  of  R according  ot  the  formula  given  in  Section  2. 

* 
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SUBROUTINE  FUNCTION  FOR  THE  UPPER  CONFIDENCE  LIMIT  OF  & = C*  - C_  • 

w 8 


100 

110 

120 

130 

140 

150 

160 

170 

130 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 


l 


FUNCTION  WU  <fi»  B»  X,  Y»  R> 

C=R 

E=B 

U=X 

V=Y 

Q=R 

QS=Q*Q 

DS=QS*  <C ♦ < 1 . +U*U> -QS> 

HL= CC^U-SORT <DS) > ✓ <C-QS> 

HU= (C^U+SQRT <DS> > s <C-QS> 

D=.  01 
X1=HL+D 

1 Y1=DU  <X1»C»E»U»V»Q) 

IF  <1 . -Yl>  2,3*3 

2 X1=X1+D 
GO  TO  1- 

3 XU=Xl-D/2. 

FR=QS-C# <U-XU> ♦ <U-XU> ✓ < 1 . +XU*XU> 
FXU=E»V/(E-FR> 

6XU-FR*E*  < 1 . +V»V> -FR^FR 
GXU=SQRT  <GXU> / <E-FR) 

FXU=FXU+6XU 

WU=FXU-XU 

RETURN 

END 

FUNCTION  DU<H,R,B,X,  Y,R> 

W=H 

C=R 

E=B 

U=X 

V=Y 

Q*R 

FR=Q»QH> (U-W> ♦ <U-W) s<l. +W»W) 

EFR=2. +C* <U-W> ♦ < 1 . +W4U) / < < 1 . +W*W> •+£) 
GFR1=E»V/ < <E-FR> ♦♦E) 

SFR=FR»E»  <1 . +V*V> -FR^FR 
GFR2=1./'SQRT<SFR) 

GFR2*  <E»  <1 . +V*V>  -2.  ♦FR)  ♦GFR2/'  <2.  ♦ CE-FR)  > 
GFR3=S0RT  <SFR> ✓ < <E-FR> ++£) 

DU=  <6FR1-M5FR2+GFR3>  *EFR 

RETURN 

END 


! 
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SUBROUTINE  FUNCTION  FOR  1HE  LOWER  CONFIDENCE  LIMIT  OF  & = £ - £ • 

t 


100 

FUNCTION  WL  <R»B»X* Y»R> 

110 

C=R 

120 

E=B 

] 

130 

U=X 

140 

V=Y 

• 

150 

Q=R 

i 

160 

QS=Q»Q 

170 

DS=QS^  < 1 . +U^U> -QS> 

180 

HL= <C*U-SGRT (DS) > ✓ <C-QS> 

■ 

190 

HU=  <C*U+SGRT (DS> ) / <C-QS> 

. 

200 

D=.  01 

1 

210 

X1=HU-D 

220 

1 

Y1=HL<X1»C*E»U.V»Q) 

230 

IF  <1 . -Yl)  2»  3»  3 

24  0 

2 

X1=X1-B 

250 

GD  TO  1 

260 

3 

XU=Xl+D/2. 

270 

FR=QS-C>  CIJ-XU)  ♦ <U-XU>  ✓ <1 . +XU»XU> 

280 

FXU=E»V'  <E-FR> 

* 

290 

GXU=FR^E»  <1 . +V^V> -FR^FR 

300 

6XU=SQRT <6XU>  ✓ <E-FR> 

i 

310 

FXU=FXU-GXU 

320 

WL=FXU-XU 

330 

RETURN 

340 

END 

350 

FUNCTION  DL<H»R,B»X-.Y*R>  ' 

360 

W=H 

370 

C-R 

380 

E=B 

j 

390 

U=X 

400 

V=Y 

\ 

410 

Q=R 

420 

FR=Q»Q-C»  CU-W) ♦ <U-W> / <1 . +W#W) 

430 

EFR=2. <U-W> ♦ <1 . +W^U) ✓ < <1 . +W^W> ♦♦E) 

440 

GFR1=E^V/'  < <E-FR>  ♦♦E) 

I 

450 

SFR=FR*E*  < 1 . +V»V> -FR^FR 

460 

GFR2= 1 . /SORT  <SFR> 

470 

GFR2=  <£♦(!.  +V»V)  -2 . ♦FR)  ♦6FR2>'  <2 . ♦ <E-FR>  ) 

480 

GFR2=-GFR2 

490 

GFR3*SQRT <SFR> ✓ < (E-FR) ♦♦E) 

500 

GFR3=-GFR3 

510 

DL= <GFR 1 +GFR2+GFR3> »EFR 

520 

RETURN 

...J 

530 

END 
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